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Abstract 



We study the first passage time (FPT) problem for biased continuous 
time random walks. Using the recently formulated framework of fractional 
Fokker-Planck equations, we obtain the Laplace transform of the FPT den- 
sity function when the bias is constant. When the bias depends linearly on 
the position, the full FPT density function is derived in terms of Hermite 
polynomials and generalized Mittag-Lemer functions. 
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I. INTRODUCTION 



Continuous time random walk (CTRW) is an important model for anomalous diffusion 
in which the mean square displacement scales with time as < X 2 (t) >~ V with < 7 < 2. 
In applications of CTRW, one is often interested in its first passage time (FPT) density 
function Q. The theory of first passage time problems has a long history PHIT] with 
applications in physics | I2]| , biology [13], and engineering [14 . 

In an earlier communication [15|, we had described the exact solution for the FPT density 
function for CTRW (in the diffusion limit) with zero bias. In this paper, we consider the 
first passage time problem for CTRW with nonzero bias (both constant and nonconstant). 
We start with a brief description of the CTRW and write down the fractional Fokker- 
Planck equation that describes it following Metzler et. al. ||16|| . Next we obtain the Laplace 
transform of FPT density function for a constant bias CTRW. Finally we obtain the full FPT 
density function when the bias depends linearly on position in terms of Hermite polynomials 
and generalized Mittag-Leflier functions [O . 



II. CTRW WITH NONZERO BIAS 

Consider a one dimensional continuous time random walk on a discrete space lattice. 
Denoting the probability of being at site j at time t by pj(t) and restricting ourselves to 
nearest neighbor jumps, the CTRW can be described by the following generalized master 



equation 16 



Pj (t) = f dt' [Aj-jpj^tf) + B j+lPj+1 (t')) 4>(t - t') + 5 Xt0 $(t). (2.1) 

Here ip(t) is the waiting time distribution, $(£) = 1 — J*dt'ip(t') is the survival probability, 
Aj-x (Bj+i) is the probability to jump from site j — 1 (j + 1) to site j. We consider a 
biased random walk where the probabilities Aj and Bj are not equal and could depend on 
the position of the random walker. However, they fulfill the local condition Aj + Bj = 1. 
The waiting time distribution is taken to be a Levy distribution |ll| with a power law tail: 



ijj(t) ~ (t/r) -1-7 for large t where < 7 < 1 and r has dimensions of time. The biased 
CTRW described above is subdiffusive i.e., its mean squared displacement varies with time 
t as t 1 where < 7 < 1. 



Following Metzler et. al. ||16|| , we can now transition to continuous space and introduce 
the Taylor series expansion: 

.. , . . . . . dA(x)p(x, t) (Ax) 2 d 2 A(x)p(x,t) 
Aj-xpj-xit) ~ A(x)p(x, t) - Ax g y + ±-J. ^ \ (2.2) 

with a similar expansion for Bj + iPj + i(t). We substitute the above expansions in the gener- 
alized master equation and Laplace transform it in time. Further, we take the generalized 
diffusion limit: Ax, r — > such that 

F(x)= lim *£[A(x) - B(x)], (2.3) 

K= i im ^1!. (2.4) 

Ax,t^o 2r 7 
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Here F(x) is the force acting on the particle during every jump and K is the generalized 
diffusion constant. Due to the presence of the force F(x), the biased CTRW is also known 
as "CTRW with (non)constant force" The Laplace transformed generalized master 

equation can be inverted (see Metzler et. al. [0 for details) to finally give the following 
fractional Fokker-Planck equation (FFPE) for the probability density function W(x, t) of 
the biased CTRW in the generalized diffusion limit: 



W{x,t) - W{x,0) 



K^W(x,t)-^F(x)W{x,t) 



< 7 < 1. 



(2.5) 



Here we have incorporated the initial condition W(x, 0) and the Riemann-Liouville fractional 



integral operator D t 7 acting on a function g(t) is defined as |T9| , |20 

1 



Drg(t) 



r( 7 ) Jo 



dt' (t-t'y^git'), 7>0. 



(2.6) 



Here T(z) is the usual gamma function pi| . We call the diffusive process obtained from 
the CTRW as "Levy type of anomalous diffusion" due to the waiting time distribution 
characteristics JT5[ . 

We now formulate the first passage time problem for the FFPE given in Eq. ( |2.5|) . This 
problem can be recast as a boundary value problem with absorbing boundaries ||. In our 
case, to obtain the FPT density function, we first need to solve Eq. ( |2.5| ) with the following 
boundary and initial conditions: 



W{0,t) = 0, W(oo,t) = 0, W(x,0) = d(x-a), 



(2.7) 



where x = a is the starting point of the CTRW, containing the initial concentration of the 
distribution. Once we solve for W(x, t), the first passage time density f(t) is given by @ 



f(t) = ~T+ / dxW(x,t). 
at Jo 



(2.8) 



III. FIRST PASSAGE TIME FOR CONSTANT BIAS 



In this section, we study the first passage time problem for Levy type anomalous sub- 
diffusion with constant bias —v (that is, F[x) = —v). The Laplace transform of the FPT 
density is obtained. 

For constant bias, the FFPE to be solved is given by [cf. Eq. ( |2.5| )] 



w(x,t)-w(x,o)= d; 



< 7 < 1. 



(3.1) 



with the boundary and initial conditions given in Eq. ( |2.7|) . We solve the FFPE using the 



method of separation of variables f22|. Let W(x, t) = X(x)T(t). Substituting in Eq. ( j3~T 
we obtain 



X(x)T(t) - X(x) = oDt 7 T(t) [uX'(x) + KX"{x)\ 



(3.2) 
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where the primes denote the derivatives with respect to x. Separating out the variables and 
introducing the separation constant A we get 

KX"{x) + vX\x) = -XX(x), (3.3) 

and 

T(t) - 1 = -A D^T(t). (3.4) 



First we solve Eq. (b\4). Taking its Laplace transform we obtain 



T(s) - - g = A T(a ). (3.5) 

Here we have used the fact that the Laplace transform of QD^Tit) is T(s)/s 7 . Solving for 
T(s) we get 

T(s) = (3.6) 

Taking the inverse Laplace transform |I7| we finally obtain 

T(t) = E, [—At 7 ] , (3.7) 
where E 7 (z) is the Mittag-Lefner function [T7]] with the following power series expansion: 



z n 



n=0 



Note that E 1 (z) reduces to the regular exponential function when 7 = 1. 

Next consider Eq. (p.3|). The solution of this equation satisfying the boundary conditions 
is given by 



sm[xJX/K -u 2 /4K 2 ] 

X(x) = exp[-v(x - a)/2K] — V . = -, A > u 2 /AK. (3.9) 

2JKX - z/ 2 /4 



Thus we have a continuous spectrum for A. Combining the solutions for X(x) and T(t), 
W(x, t) is given by the following integral over A: 



2 r°° sm\xJX/K - is 2 /4K 2 } 

W(x,t) = - / dXA(X)exp\-v(x - a)/2K] — V / -E~ [-At 7 ] . (3.10) 

kJuV±k 2 JKX - z/ 2 /4 



The coefficient A(X) is fixed by the initial condition (W(x, 0) = 5(x — a)) and we get 

2 r 00 

2 /4K 



2 r°° 1 

W(x, t) = - dX exp[-u(x - a)/2K\ sm[aJX/K - u 2 /AK 2 

IX Jv 2 /4:K V 



sm[xJX/K - v 2 /AK 2 } 

V -E 7 [-At 7 ]. (3.11) 



2 JKX - z/ 2 /4 



Letting A' = JX/K - v 2 jAK 2 we obtain 



2 r°° r i 

W(x,t) = - dX' sin X' a sin X' x ex.-p[-v( x - a) feKjE- -(KX' 2 + v 2 /AK)f . (3.12) 

Using standard trigonometric identities and dropping the primes, the above equation can 
be rewritten as 

i r°° r 

W(x, t) = - dX cos X(x - a) exp[-u(x - a)/2K\E~ -(KX 2 + v 2 /AK)t~< 

TV JO ' 1 

I C°° r i 

/ dXcosX(x + a)exp[-u(x-a)/2K]E 1 -(KX 2 + u 2 /AK)t" f . (3.13) 

7T JO 1 J 



Taking the Laplace transform with respect to time we get 

q(x,s) 



dX cosX(x — a) exp[— v(x — a)/2K] — — — r ^ — 

k Jo i s + v 2 s 1 -y/4:K + KX 2 s 1 -y 



1 r°° 

■— I dXcos X(x + a) exp[— v(x — a) 2K\ 

71 Jo 
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s + iy 2 s 1 -f/AK + KX 2 s 1 - 



(3.14) 



where q(x,s) is the Laplace transform of W(x,t). Here we have used the fact that the 
Laplace transform of E y (-BV) is l/(s + Bs 1 ^) (T?j. 

We can now perform the integration over A by using the following result [E31] 



00 cos Ax 



7T 



-a|x| 



o a 2 + X 2 2a 



(3.15) 



Thus we obtain 
1 



q(x, s) 



2y/K 



exp[— v(x — a)/2K) 



s~f + u 2 /AK 



expf-JsT + u 2 /AK \x - a\/yfK\ - exp[-Js"< + v 2 jAK (x + a)/VK] \ . (3.16) 



To obtain the Laplace transform F(s) of the FPT density function f(t), we take the 
Laplace transform of Eq. ( |2.8| ) to get 



F(s) 



dxq(x,s) + / dxW(x, 0). 



(3.17) 



Here we have used the fact that Laplace transform of dW(x, t)/dt is given by sq(x, s) — 
W(x,0). Since W(x,0) = 5(x - a) [cf. Eq. Q], we obtain 



F(s) 



1 — s dx q(x, s). 
Jo 



(3.18) 



Substituting for q(x, s) from Eq. ( |3.16| ), we get 



F(s) = 1 



dx exp[— v{x — a)/2K] 



2^fK + U 2/AK 

jexp[-yV + v 2 jAK \x - a\j\fK\ - exp[— ^/st + v 2 jAK (x + a)/V#]j . (3.19) 



After considerable manipulation, we obtain the final result upon evaluating the integrals to 
be 



F(s) 



exp 



a(y- \lv 2 + 4ATst) /2K] , < 7 < 1. 



(3.20) 



We now prove that F(s) is the Laplace transform of a valid probability distribution by 



showing that F(0) = 1 and F(s) is completely monotone p4 



; - l)n <^E!A > 0, Vn, 8 > 0. 
v ' ds n ~ 



(3.21) 



It is easily verified that F(0) = 1. We next prove that F(s) is completely monotone for 
< 7 < 1. F(s) can be written as the product Fi(s)F 2 (s) where Fi(s) = exp(ai / /2K) and 
F 2 (s) = exp(— p(s)) with 



p(s) = y/u 2 +TK&/2K. 



(3.22) 



From Feller P4| , ^(s) is completely monotone if p(s) is a positive function (for s > 0) 
with a completely monotone derivative. It is obvious that p(s) is a positive function. The 
derivative p'(s) is given by 



7s 



7-1 



V^ 2 + 4:KS^ 

Now p'( s ) can be written as the product p'i(s)p' 2 (s) where 

pf x (s) = 7s 7 - 1 , p' 2 (s) = 1/Vv 2 + 4Ks-r. 
The first factor p[(s) is completely monotone since 

(-l) n ^^ = (-l) n T(7 - 1) • • ■ (7 - n)s^- 1 > 0, 



(3.23) 



(3.24) 



(3.25) 



for < 7 < 1. Similarly, p 2 ( s ) i s a l so completely monotone. Now the product of two 
completely monotone functions is also completely monotone ||24j| . Hence p'{s) is completely 
monotone. This implies that ^(s) is completely monotone. Trivially, F\(s) is also com- 
pletely monotone. Therefore the product F(s) = F\(s)F2(s) is completely monotone for 
< 7 < 1. Consequently, F(s) is the Laplace transform of a valid probability distribution. 

From the above expressions for F(s) and its derivatives, one can calculate the moments 
of the FPT distribution [24}|. We find that all the moments diverge for < 7 < 1. 

The above solution for F(s) can be extended to include 7 = 1 (ordinary Brownian 
motion). Specifically, for 7=1, the inverse Laplace transform of F(s) can be carried out to 
give 



/(*) 



: exp 



(a - utf 
AKt 



a > 0, t > 0. 



(3.26) 



This is nothing but the FPT density function for a Brownian motion with drift v |25| . This 
can be seen as follows. For 7 = 1, the FFPE in Eq. ( p.l|) reduces to: 
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W(x,t)-W(x,0)= oD^ 1 u^- 
Differentiating both sides with respect to t we get 



r\2 

u^-W{x,t)+K—W(x,t) 



(3.27) 



-W(x, t) = »-^W(x, t) + K—W(x, t). 



(3.28) 



This is the usual Fokker-Planck equation for ordinary Brownian motion with drift v thus 
explaining the observed FPT density for 7 = 1. In other words, for a ordinary Brownian 
motion, a constant bias is equivalent to a constant drift. 



IV. FIRST PASSAGE TIME FOR LINEAR BIAS 

In this section, we study the first passage time problem for Levy type anomalous subd- 
iffusion with bias —ux which depends linearly on the position x (that is F(x) = —ujx). An 
expression for the FPT density is obtained in terms of Hermite polynomials and generalized 
Mittag-Leffler functions |17| . 

In this case, the FFPE to be solved is given by [cf. Eq. (|2.5|)1 



W(x,t)-W(x,0)= D t 



^{ U xW{x,t)) + K^W(x,t) 



< 7 < 1. 



with 



W{0,t) = W{oo,t) = 0, W(x,0) = 5(x - a). 



(4.1) 



(4.2) 



We solve this FFPE again using the method of separation of variables |[22|| . Let W(x, t) 
X(x)T{t). Substituting in Eq. ( |4.1|) we obtain 



X(x)T(t) - X(x) = \ D^T(t)} [ooxX'(x) + ooX(x) + KX"(x)\ 



(4.3) 



where the primes denote derivatives with respect to x. Separating out the variables and 
introducing the separation constant A we get 



and 



KX"{x) + tuxX'(x) + uX{x) = -XX{x) 



T{ t) - 1 = -A Q D^T(t). 



(4.4) 



(4.5) 



The above equations were solved for natural boundary conditions (iy(— oo,t) 
W(oo,t) = 0) by Metzler et. al. fll6| . We follow a similar approach here. Let 



Substituting in Eq. 



Y(x) = e~* /2 X(x), x = y/u/K x. 
and assuming to 7^ we get 



(4.6) 
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Y"(x) - xY\x) + -Y(x) = 0. (4.7) 

This is nothing but the differential equation satisfied by Hermite polynomials H n (x/y2). 
In our case we have the boundary conditions Y(0) = Y(oo) = 0. But it can be easily 
shown that any function satisfying these boundary conditions can be represented in terms 
of the odd degree Hermite polynomials H 2n +i(x/\/2). Further, these polynomials satisfy the 
following orthogonality condition which follows in a straightforward fashion from the usual 
orthogonality condition [^TJ satisfied over the domain (—00,00): 

/ dxH 2n+1 (x)H 2m+1 (x)e- x = ^2 2n {2n + l)\8 nm . (4.8) 
J 

Thus the equation for Y(x) becomes an eigenvalue equation with the eigenvalue 

A 2 n+i = (2n + IV, n = 0, 1, . . . , (4.9) 

and the eigenfunction 



Y 2n+1 {x/V2) = H 2n+1 {x/V2). (4.10) 

2 n 7r 1 /4^(2n + 1)! 



In terms of the original variables we have 



X 2n+l {x) = e-» x l 2K H 2n+l {Juj/2Kx). (4.11) 

2"vr 1 /4 J(2n+ 1)! 



We have already solved the equation for T(t) in the previous section. The only difference 
here is that A takes on discrete values given by Eq. ( |4.9| ): 

T 2n+ i(t) = E 7 [-{2n + l)ut\ (4.12) 

where Ey(z) is the Mittag-Leffler function defined earlier. The general solution for W(x,t) 
can now be written down: 

W(x,t) = ]T A 2n+1 e-^ 2/2K H 2n+1 (Juj/2Kx)E y [-(2n + 1)^]. (4.13) 

„=o 2™7r 1 /4 v /(2n + 1)! 

The coefficients A 2n+ i are determined by imposing the initial condition W(x, 0) = S(x — a). 
Using the orthogonality condition for H 2n+ i(x) given in Eq. ( |4.8| ), we can easily determine 
A 2n+ i to be 

1 



A 2n+1 = Juj/2K H 2n+1 (Ju/2Ka). (4.14) 

V 2™7r 1 /4 v /(2n + l)! V 

Substituting this in the expression for W(x,t) we get 

00 1 

W(x,t) = vWmE ^ H 2n+l {^j2Ka)e-- 2 ' 2K 



H 2n+1 Uuj/2Kx)E 1 [-{2n + l)ujf\. (4.15) 
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To obtain the FPT distribution, we substitute the above expression for W(x,t) in Eq. 
(2.8). Using the following result [21 



\ dxe-^ 2/2K H 2n+1 (Juj/2Kx) = J2K/u(-l) n ^T(n + 1/2), (4.16) 

JO v v \/7T 



we get 



fit) = -- E 1 A V J g 2n+ i(vW2^a)-E 7 [-(2n + 1)^]. (4.17) 



7T 



n=0 



r(2n + 2) z ™ v ' y dt 



From the series expansion of Mittag-Leffler function Eq. ( |3.8| ) we have 

+ 1M . | ^ | -' 2 " +i) ^f" +1) ^' 1 <-> 

Using the facts that the summand vanishes for k = and T(l + x) = xT(x) we have 

^ + - t l)"^1*^[-(2n + 1>P^1 ( , 19) 

Letting k f = k — 1 and simplifying we get 

!*,[- (an + 1 VT] = -(2n + t ^7^ - (4-20) 

But the infinite sum is the series expansion of the generalized Mittag-Leffler function 
E in [— (2n + l)^ 7 ] 01 . Therefore we have 

-^-£ 7 [-(2n + l)wt 7 ] = -(2n + l)wt 7 ~ 1 E 7i7 [-(2n + l)c^ 7 ]. (4.21) 
Substituting this in the expression for /(£) we finally get 



E i j^ r ft„ +1 (v / ^/2A'a)B 1 , T [-(2n + (4.22) 



Here we have also used the following doubling rule |2l] for gamma functions to simplify the 
expression for /(£): 

T(2n + 2) = (2n + l)r(2n + 1) = (2n + l)4 n r(n + l/2)r(n + 1)/V5r. (4.23) 



V. CONCLUSIONS 

In this paper we studied the first passage time (FPT) problem for a biased continu- 
ous time random walk (CTRW). Solution to this problem was obtained using the recently 
formulated fractional Fokker-Planck equation which describes the CTRW in the (general- 
ized) diffusion limit. Expressions for the FPT density function was obtained in two cases - 
constant bias and bias that depends linearly on the position. 
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